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Abstract 

We discuss the first applications of our newly developed Monte Carlo event gen¬ 
erator SOPHIA to multiparticle photoproduction of relativistic protons with thermal 
and power law radiation fields. The measured total cross section is reproduced in 
terms of excitation and decay of baryon resonances, direct pion production, diffractive 
scattering, and non-diffractive multiparticle production. Non-diffractive multiparti¬ 
cle production is described using a string fragmentation model. We demonstrate that 
the widely used ‘A-approximation’ for the photoproduction cross section is reasonable 
only for a restricted set of astrophysical applications. The relevance of this result for 
cosmic ray propagation through the microwave background and hadronic models of 
active galactic nuclei and gamma-ray bursts is briefly discussed. 

Keywords: nuclear reactions - elementary particles - radiation mechanisms: non- 
thermal 


1 Introduction 

The energy spectrum of cosmic rays extends to energies above 10 20 eV. Such particles are most 
likely extragalactic and propagate through many tens of Mpc before reaching the Earth. The 
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main energy loss mechanism for nucleons in this energy range is photomeson production on 
the cosmic microwave background radiation (CMBR). This process has an energy threshold 
of 1.08 GeV in the center of momentum frame (CMF) of the interacting particles. It causes 
the distortion of the proton spectrum above 3 x 10 19 eV during propagation, known as the 
Greisen-Zatsepin-Kuzmin cutoff (Greisen 1966, Zatsepin & Kuzmin 1966; see Protheroe & 
Johnson 1996 for additional references). In environments harbouring dense radiation fields 
with higher photon energies the particle production threshold is reached at lower proton 
energies, e.g. ~ 10 16 eV for ambient photon energies of 10 eV. If protons are accelerated in 
such energetic astrophysical objects, photopion production and subsequent pion decay would 
lead to the emission of energetic y-rays and neutrinos. This may be observable from jets 
in active galactic nuclei (AGN) (Mannheim 1993; Protheroe 1996) and gamma-ray bursts 
(GRB) (Waxmann & Bahcall 1997; Vietri 1998a,b; I Pit teller & Dermer 1998; Rachen & 
Meszaros 1998). 

The prominent A(1232)-resonance near the threshold has often been used to construct 
approximate pion production cross sections (referred to as the ‘A-approximation’ hereafter) 
and to determine gross features like y-ray-to-neutrino energy yields, proton inelasticities, etc. 
which enter the relevant astrophysical calculations. In this paper we discuss first results using 
a new Monte-Carlo event generator for photohadronic interactions of relativistic protons in 
radiation fields of astrophysical origin, which includes all relevant interaction processes and 
is based on models and data available from particle physics. 

Our event generator, SOPHIA (Simulation Of Photo-Hadronic Interactions in Astro¬ 
physics), has been extensively tested on all fixed target and collider experiment data avail¬ 
able to us (Miicke et al. 1998a). A comparison of the results from SOPHIA with the A- 
approximation emphasizes the importance of including all interaction processes in calcula¬ 
tions of expected astrophysical signals. Detailed estimates of such signals will be subject of 
future work. This paper concentrates on photopion production features that may be impor¬ 
tant for astrophysical processes which involve hadronic reactions. In Sec. 2 we summarize the 
physics implemented in SOPHIA and contrast it with the A-resonance approximation. Var¬ 
ious astrophysical applications of SOPHIA are discussed in Sec. 3 by considering power-law 
and thermal photon seed spectra. A summary is given in Sec. 4. 

2 Cross section and kinematics 

Consider a relativistic proton with energy E p = y m p and rest mass m p , which interacts 
with a photon of energy e at an angle 6. The square of the total CM frame energy of the 
interaction is given by 

s = rrip + 2E p e(l — (3 cos 0) = vr? v + 2 m p e', (1) 

where j3 = a/1 — y -2 is the velocity of the proton in terms of the velocity of light and 
e' the photon energy in the nucleon rest frame (NRF). In SOPHIA, proton and neutron 
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induced interactions are distinguished - for simplicity we quote only proton interactions in 
the following. 

The partial cross sections for resonance excitation, direct (non-resonant) single-pion pro¬ 
duction, and diffractive scattering are determined by fits to exclusive data (Fig. 1). We 
consider the 9 most important resonances (A + (1232), iV + (1440), iV + (1520), iV + (1535), 
iV + (1650), iV + (1680), A + (1700), A + (1905) and A + (1950)) and use the Breit-Wigner for¬ 
mula together with their known properties of mass, width and decay branching ratios to 
determine their contribution to the individual interaction channels. The direct channel is 
defined as the residual, non-resonant contribution to the channels py —> mi + , py —> A ++ n ~, 
and py —> A°7 t + . The method is described in detail in Rachen (1996). Note that, although 
the A(1232)-resonance has the highest cross section at low interaction energies, the direct 
channel dominates near threshold. The direct channel can be important for proton interac¬ 
tions in soft photon spectra since it produces exclusively charged pions. Diffractive scattering 
is due to the coupling of the photon to the vector mesons p° and u, which are produced at 
very high interaction energies with the ratio 9:1 and a cross section proportional to the 
total cross section. Finally, the residuals to the total cross section are fitted by a simple 
model and treated as statistical multipion production which is simulated by a QCD string 
fragmentation model (Andersson et al. 1983). After the decay of all intermediate states, 
considering basic kinematical relations and accelerator data on rapidity distributions, the re¬ 
sulting distributions of protons, neutrons and pions in the source determine the measurable 
astrophysical quantities, for example, the cosmic ray, neutrino and y-ray emission. Details 
of the simulation techniques are described in Miicke et al. (1998a). 

In the A-approximation (Stecker 1973, Gaisser et al. 1995) the total cross section is given 
as ua = 500 pbarn ©(a/s — niA + r A /2) • @(m A + Ta/2 — y^), where = 1.232 GeV is the 
mass and Ta = 0.115 GeV is the width of the A(1232)-resonance, and 0 is the Heaviside 
step function. The A-approximation uses the branching ratios of the A + (1232)-resonance 
to determine the number ratio ir° to n + of 2:1. Photohadronic z/-production is the result 
of the decay of charged secondary pions (7r + —> e + v p v p v e , 7—> e _ z/ /t z> M P e ). Gamma rays 
are produced via neutral pion decay ( 7 r° —» yy) and synchrotron/Compton emission from 
the resulting relativistic leptons. In the A-approximation this leads to a ratio of the energy 
content in y-rays to neutrinos E 1 : E u = E 1 : E u — 3:1. The decay kinematics of this 
resonance decay predicts a nucleon inelasticity K p = A E p /E p & 0.2. 

3 Astrophysical applications 

In this section we illustrate the importance of using the complete photohadronic cross section 
and final state composition and kinematics for a number of astrophysical applications. In 
particular, we compare the average proton inelasticity in a photopion production K p and the 
ratio E 1 jE v from simulations with SOPHIA to the A-approximation. We assume that due 
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Figure 1: The total /ry cross section with the contributions of the baryon resonances consid¬ 
ered in this work, the direct single-pion production, diffractive scattering, and the multipion 
production as a function of the photon’s NRF energy (1 //barn = 10 -30 cm 2 ). Data are from 
Baldini et al. (1988). 


to radiative processes electrons convert all their kinetic energy into photons, and so e ± are 
counted in £ 7 . The total neutrino energy £ u is the sum of the energies of z/ e , v p and their 
antiparticles. 

3.1 Power—law seed photon spectra 

In order to illustrate the contribution of the different interaction processes in various astro- 
physical environments we convolute the cross section with proton and seed photon spectra. 
Proton spectra as well as non-thermal seed photon spectra can be well approximated by 
power laws. Figure 2 shows the distribution of the number of interactions with CMF energy 
squared, s, for power-law seed photon spectra n e ~ e~ ai and a proton spectrum n® ~ E~ ap 
with index a p = 2. Spectra with a p = 2 are typical for shock accelerated particles. Shock 
acceleration is believed to be the dominant process for acceleration of galactic cosmic rays 
(Jones & Ellison 1991). It may also be present in extragalactic particle accelerator (Bicr- 
mann & Strittmatter 1987), in particular AGN (e.g. Stecker et al. 1991, Mannheim 1993, 
Protheroe 1996) and in GRBs (e.g. Waxman 1995, Vietri 1998). 
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Figure 2: Interaction probability distribution as function of the squared CMF energy with 
index a 7 and proton spectrum ng ~ E~' 2 . 

It has been proposed that the high energy gamma ray emission from blazar jets as 
observed in the GeV to TeV energy range results from photopion production in a proton- 
electron plasma where the relativistic electrons provide a synchrotron seed photon held for 
the high energy protons (“Proton Blazar” model, Mannheim 1993). The synchrotron photon 
spectrum is described by a power law which is often observed in jets of AGN. Photopion 
production is followed by electromagnetic cascades reprocessing the injected power of the pion 
decay products (Mannheim et ah 1991). Therefore, Fig. 2 may approximate the situation in 
AGN jets. While for steep photon spectra the prominent resonance region near threshold is 
the dominating contributor to the interaction rate, the high energy region of the cross section 
gains importance for flat seed photon spectra. Radio spectra with a 7 « 1 (F v oc const.) up 
to a break frequency v b are, for example, observed in flat spectrum radio quasars (FSRQs) 
and BL Lac objects. Such spectra are believed to originate from a superposition of several 
self-absorbed synchrotron components (Cotton et al. 1980, Shaffer & Marscher 1979). Above 
u b the spectrum is loss dominated (i.e., the cooling time scale is shorter than the dynamical 
time scale) and steepens to a 7 ~ 2. Recent observations of blazars have revealed synchrotron 
peaks ranging from the optical/IR-band up to UV/soft X-rays, especially in the flaring 
state (see, e.g., Fossati et al. 1998). Thus flat power law photon spectra exist in AGN 
jets and would shift the average CMF energy of photopion production towards high values. 
For protons with energies below 10 16 eV only the infrared and higher frequency photons 
are relevant for photopion production due to the threshold condition. For blazars with 
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a synchrotron peak at rather low frequencies the target photon spectrum for photomeson 
production may be rather steep, and photohadronic interactions would mainly occur near 
the threshold. 

In both, flat and steep photon spectra, the simulations with SOPHIA show that the 7 - 
ray-to-//-energy ratio is £ 7 /£ v ~ 1. The difference with the A-approximation {S 1 /S v = 3) is 
due to the contribution of the second resonance and the multipion production regions in the 
case of flat photon spectra. In steep spectra, the difference is explained by the influence of 
the immediate threshold region, where the direct channel dominates over the Delta resonance 
(Mhcke et al. 1998b). Therefore, the 7 -ray and neutrino outputs are approximately equal. 
This may have significant impact on estimates of the predicted neutrino fluxes from AGN 
jets which are usually normalized to the observed 7 -ray luminosity. 

The absolute 7 -ray and neutrino energy outputs depend on the total proton energy loss 
K p E p . For photon power law indices a 1 > 2 we find K p « 0.2, as in the A-approximation. 
The inelasticity increases rapidly for flatter photon spectra to reach K p « 0.6 for a 1 = 1. 
This shows that the A-approximation underestimates significantly the fractional energy loss 
of the incident nucleon in environments with flat seed photon spectra, where protons cool 
much faster. 

Proton energy losses, however, limit the maximum possible proton energy which can be 
achieved in AGN jets (Biermann & Strittmatter 1987). Setting the time scales for proton 
acceleration and energy loss equal, Mannheim (1993) obtained cutoff energies of ~ 10 16 eV. 
This estimate is based on a mean cross section for photopion production ( K p a pl ) = 50/ibarn 
and a mean inelasticity (K p ) = 0.25. Noting that for flat seed photon spectra SOPHIA 
simulations give proton inelasticities about a factor of 2 larger than Mannheim (1993), we 
expect slightly lower maximum proton energies. Similar corrections could apply to any 
proton acceleration model where the maximum proton energy is limited by photohadronic 
interactions. In contrast, models in which Larmor radius constraints and finite acceleration 
time scales limit the proton energies, would not be affected (e.g. proton acceleration in hot 
spots of radio galaxies, Rachen & Biermann 1993). 

An ideal environment for proton acceleration may be found in gamma-ray bursts which 
have a broken power-law photon spectrum steepening from a 7 « 1 to a 7 ~ 2 at keV 
energies (in the comoving frame). Some GRB models involve the acceleration of nltrahigh 
energy cosmic rays. Cascades initiated by photopion production follow (Waxman & Bahcall 
1997, Vietri 1998, Bottcher & Dermer 1998, Rachen & Meszaros 1998). Figure 3 shows 
the probability of photoproduction interaction at squared CMF energy s for a iW 2 -proton 
spectrum with a photon spectrum with a 7 = 2/3 below 1 keV and a 7 = 2 above that 
energy. The curves are calculated for different parts of the proton spectrum, extending from 
minimum energies of 10 12 —10 18 eV to a high energy cutoff of 10 21 eV. While for a proton 
spectrum with E m in = 10 12 eV most interactions occur on the e~ 2 -part of the photon held, 
for higher E m \ n proton-photon interactions involve the hat part of the seed photon spectrum 
leading to higher average s values. Consequently, a 7 -ray-to neutrino energy ratio of 1 and 
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Figure 3: Interaction probability distribution as function of the squared CMF energy for a 
broken power-law seed photon spectra as observed in GRBs (n e r^j £ 2 /3 f or 10 8 eV < e < 1 
keV and n e ~ e -2 for 1 keV< e <100 keV and proton spectrum tie E 2 with E m \ n < E < 
10 21 eV and different minimum energies E m ; n = 10 12 eV (solid line), 10 14 eV (dashed line), 
10 16 eV (dashed-dotted line) and 10 18 eV (dashed-dashed line). 


inelasticities of the order of 0.5—0.7 are expected for the production of the highest neutrino 
and secondary photon energies (Miicke et ah, 1998b). 

Recently, Vietri (1998a,b) has claimed that in situ photoproduced neutrinos in the blast- 
waves associated with GRBs can reach energies up to «10 19 eV. In order to calculate the 
maximum neutrino energy, he assumed that individual neutrinos typically carry 5% of the 
incident proton energy which is derived from the A-approximation. However, SOPHIA sim¬ 
ulations show that this value is rather about 1% due to the increasing neutrino multiplicity, 
constraining the maximum neutrino energies to lower values (Miicke et al. 1998b). 

Since the average interaction energies are well beyond the resonance region in GRB we 
expect that detailed calculations with SOPHIA will change significantly the predictions of 
hadronic GRB models. 
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3.2 Thermal seed photon spectra 

By contrast to power-law photon spectra, thermal (black body) photon fields are strongly 
concentrated about a characteristic energy kT. This concentration of the seed photons leads 
to an emphasis of one specific photoproduction energy range. We consider the interactions 
of protons with different power law spectra with photons from T = 2.73 K and kT = 10 eV 
blackbody spectra. The former case is of particular astrophysical interest for cosmic ray (CR) 
propagation through the cosmic microwave background. The latter case may approximately 
describe the situation in AGN jets with the accretion disk photons being the target for 
proton photon interactions. We show the interaction probability as a function of s in Fig. 4. 
In the case of CR transport through the CMBR the A(1232)-resonance region plays the 
dominant role. SOPHIA simulations give a photon-to-neutrino energy ratio £ 7 /£7 ~ 1.7 
and inelasticities ~ 0.2. We note that the neutrino component dominates directly at the 
threshold with £ 1 /£ u ~ 1/2 due to the dominance of the direct channel. 

Recent AGASA measurements seem to indicate an extension of the CR spectrum beyond 
10 20 eV (Takeda et al. 1998). It is generally assumed that acceleration scenarios are limited 
to maximum energies of 10 21 eV. If the CR proton spectrum extends above 10 21 eV, one may 
have to consider so called top-down (TD) models where the highest energy cosmic rays are 
decay products of GUT scale ultraheavy (10 25 eV) particles (e.g., Hill 1983; Sigl et al. 1994; 
Protheroe & Stanev 1996). While for TD models the photon-to-neutrino energy ratio would 
reach unity and proton inelasticities can rise up to 0.6, in case of acceleration scenarios the 
A-approximation describes sufficiently well the proton propagation through the CMBR. 

In previous Monte Carlo calculations of cosmic ray pion photoproduction interactions 
(e.g. Protheroe & Johnson 1996), an earlier photoproduction event generator based on the 
observed inclusive cross sections for py —> n°p, and for py —> 7r + n, together with a semi- 
empirical model for multi-pion production (with equal numbers of 7r + , n~ and 7r° produced 
in the central region, and charge ratios in the fragmentation regions based on naive quark 
model arguments; Szabo & Protheroe 1994) was used. The results from this previous event 
generator are not too different, and give £ 7 /£„ ~ 1.6 and inelasticity «0.2 at the peak of the 
resonance, and £ 1 /£ u ~ 1.0 and inelasticity ~0.55 in the extreme multipion region. 

For photopion production in an isotropic UV radiation field, the high energy and the 
resonance regions of the cross section are of equal importance (see Fig. 4). Thermal radiation 
peaking in the UV-range of the electromagnetic spectrum is observed from some radio- 
loud AGN (for example 3C273). Therefore a hadronic blazar model has been proposed in 
which protons accelerated in some regions of the jet interact with accretion disk photons 
(“Accretion Disk Proton Blazar Model”, Protheroe 1997, 1998). In this model the protons, 
being charged, are isotropised in the jet frame and readily interact with the almost black-body 
photons coming from the inner part of the accretion disk near the base of the jet, provided 
the disk is sufficiently luminous or the emission region is not too far away from the disk. 
The model is tuned such that the resulting y-rays, which are relativistically beamed along 



the jet direction, undergo a pair-synchrotron cascade in the accretion disk radiation and jet 
magnetic held such that the emerging radiation can reproduce features of the observed 7 -ray 
emission of blazars. The physics implemented in SOPHIA allows for a detailed investigation 
of the parameter space of this model. 



Figure 4: Interaction probability distribution as function of the squared CMF energy for 
black body seed photon spectra and proton spectra with index a p between 10 6 , " 12 GeV. 


4 Conclusions 

We use our recently developed Monte Carlo code for photohadronic processes (SOPHIA) 
to compare predictions of astrophysically relevant quantities with the expectations from 
the widely used ‘A-approximation’. We discuss proton inelasticities and 7 -ray-to-neutrino 
energy ratios in the framework of astrophysical environments like AGN, GRBs and for the 
case of CR propagation. SOPHIA uses the full information of the cross section and the 
final state particle composition and kinematics of the interaction processes as provided by 
particle physics. We show that the A-approximation is a reasonable approximation only for 
a restricted set of applications. Several astrophysical cases are discussed: 

• Photohadronic interactions of high energy nucleons in typical radiation fields of GRBs, 
which may be responsible for the highest energy 7 -rays and neutrinos, occur mainly 
at high CMF energies. Nucleon inelasticities of « 0.5—0.6 are expected in this case 
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whereas the A-approximation gives only ^0.2. The 7 -ray and neutrino energy outputs 
are approximately equal. 

• High-energy 7 -ray emission from blazar jets may be explainable by photopion produc¬ 
tion of relativistic protons bathed in a synchrotron or thermal (UV-photons from the 
accretion disk, IR-photons from the molecular torus) radiation held (“Proton Blazar” 
models). For a hat-spectrum seed photon held photohadronic interactions with high 
CMF energies become important. SOPHIA simulations show that protons may cool 
much faster than in the A-approximation because of the growth of the inelasticity K p . 
The 7 -ray-to-neutrino energy ratio may approach unity for both, hat and steep seed 
photon spectra, whereas many models use E 1 IE V = 3. This would affect estimates of 
the predicted neutrino hux. 

• Interactions of ultra-high energy cosmic rays propagating through the microwave back¬ 
ground mainly happen in the region of the A(1232)-resonance for energies <10 21 eV. 
Above that energy the effects of the second resonant and multiparticle regions increase 
the fractional energy loss and alter the hnal shape of the cosmic ray spectrum. 

Detailed calculations addressing photomeson production in the above mentioned astrophys- 
ical environments will be reported in forthcoming publications. 
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